Tests and applications of self-consistent cranking in the interacting boson model 
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The self-consistent cranking method is tested by comparing the cranking calculations in the inter- 
acting boson model with the exact results obtained from the SU(3) and 0(6) dynamical symmetries 
and from numerical diagonalization. The method is used to study the spin dependence of shape 
variables in the sd and sdg boson models. When realistic sets of parameters are used, both models 
lead to similar results: axial shape is retained with increasing cranking frequency while fluctuations 
Q"^ , in the shape variable 7 are slightly reduced. 
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I. INTRODUCTION 



Self-consistent cranking (SCC) is one of the most popular methods to study collective rotations in nuclear many- 
body systems [jjj. Within the Hartree-Fock framework, it has provided important insights into the backbending 
phenomena, and currently, it is being used actively in studies of superdeformed nuclei In stark contrast, it has 
been rarely used in the interacting boson model (IBM) Q. In fact, there is only one application of SCC to IBM where a 
study of moment of inertia (MOI) is carried out Q . In the early stages of the model, an obvious reason for this neglect 
was the availability of exact results (either via dynamical symmetries or numerical diagonalization), which left little 
room for development of approximate methods such as SCC. In later extensions of the IBM (e.g. sdg-lBM), dynamical 
Q\ , symmetries were not very useful, and numerical diagonalization could not be carried out due to large basis space, 
hence alternative methods were needed. At around the same time, however, an exact angular momentum projection 
technique was developed for axially symmetric boson systems (1/N expansion |||), which provided analytical solutions 
for general Hamiltonians in large basis spaces, filling the gap left by the other methods. Although the assumption 
of axial symmetry is reasonable for most deformed nuclei in their ground states, triaxial effects are known to play 
a role, especially at high-spins Q. Therefore a study of the evolution of shapes in the IBM would be useful in 
order to investigate such questions as the effects of triaxiality in nuclear spectra, ways of incorporating it in the IBM 
Hamiltonian, and whether an improved description of spectra can be obtained if triaxial effects are taken into account. 

In this paper, we first present tests of the SCC by comparing the cranking calculations with the exact results 
obtained in the SU(3) and 0(6) limits, as well as from diagonalization of more realistic Hamiltonians in the sd-IBM. 
The SCC method is then used in a study of triaxiality in the IBM both in the sd and sdg-boson versions of the model. 



II. CRANKING FORMALISM IN THE IBM 

The SCC in the IBM was formulated in Ref. [0] to which we refer for details, especially concerning the construction 
of the intrinsic state and the symmetries it possesses. For a given IBM Hamiltonian H and intrinsic state \N, x), the 
SCC around the x-axis is described by 

8{N,x\H-u>L x \N,x) = 0, (1) 

where L x is the x-component of the angular momentum operator, N is the boson number and x are the variational 
parameters to be determined from the extremum condition. For convenience, we consider a general formulation of the 
IBM which will be tailored to specific cases later. Thus, we introduce the boson creation and annihilation operators 
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b\ , with I = 0, 2, 4, . . . where bo — s, bi — d.b^ — g, etc. The intrinsic state is given by a condensate of intrinsic 
bosons as 

|JV,x) = (N\)-^(tf) N \0), ftt = J2 £ x lm b\ m , (2) 

/ TCI — — l 

where xi m are the normalized boson mean fields, i.e. x ■ x = 1. Due to the symmetries in the cranked system [Ej, 
xi m are real and x;_ m = xi m . Note that the odd-m components of Xi m vanish in the static limit (lo — 0), but not 
in general. Thus, invoking the normalization condition (e.g. setting xqo — 1), there are 3 independent variational 
parameters in the sd-IBM (X20, %n, X22) and 5 more in the sdg-lBM (2:40, . . . , £44). For the IBM Hamiltonian, we use 
the generic multipolc form 

# = K fc T«.T« (3) 

Z fc=0 

where the boson number and multipole operators are given by 

^ = E 6 lA- rW = X)t fcil [ftt6 I ]W. (4) 

In particular, the angular momentum operators are 

L X = ^(L +1 -L_ X ), L p =X) V '^ + l)(2i + l)/3 (5) 

The parameters of the Hamiltonian consists of the single boson energies, the multipole coupling strengths, Kk, and 
the multipole parameters, tkji- The hermiticity condition on multipole operators requires that tkji — tkij- 

The expectation values in Eq. @ can be evaluated in a straightforward manner using boson calculus techniques. 
For the one-body operators, one simply obtains 

(N,x\hi\N,x)=N^2^. (6) 

m 

Calculations for the multipole interactions are somewhat more involved but the expectation values can again be 
reduced to a compact form 



(iV,x|TW.T( fe )|7V,x) = iV(iV-l) {-lYA^Ak-^ + N^'i^- (?) 

fj,= — h Im 



Here corresponds to the expectation value of the operator in a single boson state and is given by 

Ak, = £ (-ir{jnlm\kn)t kjl ^^. (8) 

jnlm 

Using the symmetries imposed on the mean fields and the hermiticity condition, it can be shown that Ak^ vanishes 
if k + /i is odd, and furthermore = (— 1) Ak^- Thus, for even multipoles, only the positive, even /i values 

contribute to the sum in Eq. (0). The second term in Eq. (Q) arises from the normal ordering of the boson operators 
and corresponds to an effective one-body term with boson energies 

, _ 2k + 1 



T Eif ( fl ) 

Finally, the expectation value of L x is 
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(N,x\L x \N,x) = -Nj2(-^) m (lnl - m\ll) y/2l(l + 1)(2Z + l)/3 (10) 
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For a given angular frequency u>, the cranking equation (|l|) is solved numerically by minimizing the cranking 
expression E(x, oj) = (H — loL x ) with respect to the mean fields x. The dynamical MOI is calculated from the 
derivative of (L x ), Eq. @ as @ 



r(2) _ d(L x 



(11) 



which is obtained under the assumption (L x ) ~ L. For the exact energy levels E(L), dynamical MOI and the 
corresponding rotational frequencies are calculated from H 



r d 2 E(L) 
dL 2 



Ah 2 



E(L + 2) - 2E(L) + E(L - 2) ' 



u = —[E(L + 2)-E(L-2)]. 



(12) 



Evolution of triaxiality with rotation is studied using the following expressions for the collective shape variable 7 



tan 7 = V2 y— 7, 
wo/ 

„ [f (Q-Q-Q) 



(13) 
(14) 



The one-body expectation values of the quadrupole operator in Eq. (|T^) follow from Eq. (||) as (Q/j) = NA^. The 



two-body expectation value in Eq. ( |l4| ) is given in Eq. (|7) and the three-body part is derived in Appendix A. Eq. (|l_ 
gives the average value for 7 while Eq. ( Jl4| ) probes its fluctuations from this average. Thus the information content 
of the two expressions for 7 are different and they compliment each other 



III. CRANKING IN THE sd-IBM 



In this section, we perform cranking calculations for a variety of sg?-IBM Hamiltonians and compare the results 
for the dynamical MOI and 7 with the exact ones obtained from dynamical symmetries (SU(3) and 0(6)), and from 
numerical diagonalization. For this purpose, we use a simple Hamiltonian with a quadrupole interaction and d-boson 
energy 



H = —kQ ■ Q + eh c i- 



(15) 



Following the convention, we denote the parameters of the quadrupole operator by £202 = 1 and £222 = X- The SU(3) 
limit is obtained when e = and x = — V< /2, and the 0(6) limit when e = x = 0- To test more realistic Hamiltonians, 
we use the middle value for \ (~ V7 /4), which leads to a better description of electromagnetic properties Inclusion 
of the ci-boson energy with e = 1.5Nk improves, in addition, the moment of inertia systematics ||. The exact energy 
eigenvalues E(L) are obtained from the Casimir operators in the case of the dynamical symmetries, and by numerical 
diagonalization of the Hamiltonian in the latter two cases. 

For convenience in the variational problem, we set xoo = 1 and drop the subscript 2 from the quadrupole mean 
fields, i.e. X2 m = x m . The normalization is then given by 



N = 1 + x 2 + 2x\ + 2irJ- 

The expectation value of the Hamiltonian jl5|), E(x) = (N,x\H\N,x), can be written from Eqs. as 

E(x) = k [N(N - 1) (A 2 Q + 2A 2 2 ) + 5(n s ) + ( X 2 + l)(n d )] + e(n d ). 
Here the quadratic forms follow from Eq. (g) as 

1 



(16) 
(17) 
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2x - ^2jl X (4 +x 2 - 2x 2 ) 

2x 2 + Vv7x (2x0x2 + Vyzxi 



(18) 
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and the expectation values of fi s and fid are 

f* \ N 

From Eq. ( |l0| ) , the cranking term is given by 



\ n d) 



TV 



2xi 



2x\) 



(L x ) = (VQxq + 2x 2 



(19) 



(20) 



In Fig. |, we compare the cranking results for the dynamical MOI obtained from Eq. ([fl]) with those obtained from 
Eq. ( fl2| ) using the exact energy levels for TV — 10 bosons. In the SU(3) and 0(6) limits, the MOI is constant and 

(2) 2 

given by Jcx = 4/3k = 200/3 and 2/k — 100 fi /MeV, respectively Note that both results are independent of TV, 
and the constancy of MOI simply follows from the fact that E{L) is quadratic in L in both cases. The SCC leads to 
a constant MOI in the SU(3) limit, as expected from a rigid rotor, but deviates from the exact result by about 5%. 
The 0(6) limit corresponds to a 7-unstable rotor for which cranking is not expected to work well. Not surprisingly, 
comparison of the cranking results for J^> in Fig. [l] indicate a maximum deviation (~ 10%) from the exact result, 
as well as a small frequency dependence in MOI. We will comment on these deviations further when discussing the 
TV dependence of the results below. The case with x — — V7/4 falls in between the two dynamical symmetry results. 
The dynamical MOI increases with w as in a typical deformed nucleus (Fig. Q), though the overall magnitude is still 
too large. The deviation between the cranking and the exact results in this case is comparable to that in the SU(3) 
limit, with the error getting smaller at higher frequencies. When a <i-boson energy with e = 1.5Nk is included in 
the Hamiltonian, which fits the observed range of MOI data better, the agreement between the cranking and exact 
results improve markedly, especially at lower frequencies. 

The SCC is a semi-classical theory and it would be exact in the classical limit when TV — ► 00. Thus to understand 
the nature of discrepancies seen in Fig. [[] better, we need to study the TV dependence of the results. In Fig. ||, we 
show the SCC results for as a function of TV in the SU(3) and 0(6) limits. Since the MOI is not constant in 
the 0(6) limit, we have taken the value of at the middle-frequency for each TV. The curves that trace the SCC 
results are obtained from 



in the SU(3) limit, and from 



J, 



(2) 



J, 



(2) 
ex 
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2TV 



7TV 2 



(21) 



7(2) = 7(2) 



TV 2 TV 2 



(22) 



in the 0(6) limit. The coefficients in E 
(0). It is clear from Fig. | and Eqs. (§1 



js. ( |2l| ) and ( p2[ ) are derived from a 1/TV expansion of the SCC equation 
f22|) that the SCC calculation of J^ 1 is correct to leading order but fails in 
higher orders in 1/TV by generating spurious correction terms. A comparison of Eq. (|l7]) with the accurate angular- 
momentum-projected result II confirms that the incorrect treatment of the 1/TV terms in the SCC is the cause of the 
discrepancy. It is interesting to note that an agreement is obtained if one uses the Casimir operator of the SU(3) given 
by C2 = Q- Q + (3/8)L • L, instead of just Q ■ Q. In this case, the incorrect 1/TV contribution from the Q ■ Q interaction 
is exactly canceled by the spurious contribution from the L ■ L term. Another interesting observation is that the best 
agreement of SCC with the exact results is achieved in the realistic case with \ = ^^7/4 and e — 1.5TV/C. This 
happens because the 1/TV term is now dominated by the leading contribution from {fid), which is correctly treated in 
the SCC. The errors from (fid) contribute at the 1/TV 2 level, and these become apparent only at the high-rotational 
frequencies as is apparent in Fig. |l|. 

We next study the effects of triaxiality in SCC. It is well known that the energy surface of an sd-lBM Hamiltonian 
with one- and two-body interactions has an axially symmetric minima in the deformed phase [fl0|| . An exception 
occurs in the 0(6) limit, where a 7-unstable shape develops. The SCC calculations in the 0(6) limit give 30° for both 
the average 7 ([jj]) and its fluctuations (|l4|) at all frequencies. Thus the SCC results are consistent with the description 
of the 0(6) limit using a triaxial intrinsic state with 7 = 30°, which has a very shallow minimum in the 7 direction 
p| . Frequency independence of the results corroborates with the fact that the 0(6) states remain 7-unstable at all 
spins. In Fig. |3[ we show the evolution of the triaxiality angle 7 with the cranking frequency for the remaining three 
cases discussed above for three different boson numbers, N — 10, 15, 20. Results obtained from both Eq. ( |l3| ) (solid 
line) and (|l4|) (dashed line) are shown. The ground band in the SU(3) limit can be described exactly by angular 
momentum projection from an axially symmetric intrinsic state, therefore, both average 7 and its fluctuations should 
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vanish at all frequencies. The SCC results (Fig. left) start with 7 = at low u but deviate from it systematically 
with increasing frequency. Although the situation appears to be improving with increasing N, similar to the case 
in J^ 2 \ in fact, this is merely due to the extension of the spectrum with N (w max oc N). If one scales out the N 
dependence (i.e., plots 7 against uj/N), then all the curves with different N overlap. Thus 7 can be written as a power 
expansion in lo/N with the leading (N independent) term being zero. This situation is similar to that encountered 
in the study of J^ 1 : the SCC gets the leading order term correctly (7 = 0) but fails in higher-orders in 1/N by 
generating spurious terms in powers of lu/N. The fact that the incorrect 1/N terms in 7 all depend on the cranking 
frequency makes the SCC results increasingly unreliable at higher frequencies. The onset of triaxiality observed at 
around mid-frequency (w max /2) in all curves in Fig. || is thus due the failure of the SCC and not a genuine feature of 
the boson system. We remark that despite the sudden increase in the calculated 7 values, remains constant in 
the SU(3) limit, suggesting that triaxiality has a negligible effect on MOI. 

We interpret the remaining more realistic cases in the light of the SU(3) test. In the case of \ = ~V7 /4 (Fig. [|, 
middle), the average 7 (solid lines) remains around zero except at high frequencies. Attributing this sudden rise in 
7 to the break down of SCC, we see that 7^0 for all ground-band states of a Q ■ Q Hamiltonian. The fluctuations 
(dashed lines), on the other hand, are non-zero but get smaller with increasing N. The values of fluctuations at u> = 
are consistent with those obtained from an exact calculation for the ground state ]T^| . At higher frequencies there 
is a small reduction in fluctuations, that is, rotations have a slightly stabilizing effect on the 7 motion. Increasing 
the boson number also suppresses fluctuations, which is due to the energy surface becoming deeper in the 7 direction 
with larger N values. The last case with e = l.hNn (Fig. ||, right) has broadly the same features. The one-body term 
cannot induce triaxiality alone, therefore, the small deviation in average 7 from zero at mid-frequencies is presumably 
due to errors in SCC. There is an increase in fluctuations compared to the middle panel, which can be explained as 
due to the energy surface becoming shallower with the addition of the one-body term. These results indicate that 
rotations have a rather limited effect on the shape variables in the sd-IBM, with average 7 remaining nearly constant 
at zero, and its fluctuations being slightly reduced from the ground state value at higher spins. 



IV. CRANKING IN THE sdg-TBM. 



In contrast to the sd-IBM, where numerical diagonalization is a routine task, the sdg-lBM already suffers from the 
large basis problem, and exact diagonalization is not possible for N > 11, that is, for most of the deformed nuclei. 
The 1/N expansion circumvents this problem but because it assumes axial symmetry, one cannot use it to address 
questions on triaxiality. Here we perform cranking calculations in the sdg-IBWL to study the evolution of shapes with 
rotation without restriction to axial symmetry. For the Hamiltonian, we choose 



H = -k 2 Q ■ Q - kaT a ■ T 4 + e d h d + e g h g , 

which has been shown to provide a good representation of data in rare-earth and actinide nuclei JlJ 
III, we set 502 — h-02 = 1 and xqq — 1. Then the normalization is given by 
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The expectation value of the Hamiltonian (23) follows from Eqs. ( 
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2x 



43 



2x 



44- 



(23) 

As in Section 
(24) 



E(x) = -K2 [N(N - l)(Aj + 2A 2 22 ) + 5{n s ) + (1 + q\ 2 + q 2 24 )(h d ) + (5/9)(g| 4 + q 2 44 )(h g )] 
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N(N - l)(Al + 2A\ 2 + 2A 2 4i ) + 9(h s ) + (9/5)(h 2 22 + h 2 24 )(n d ) + (1 + h% A + hl 4 )(n g ) 



+£d(n d ) + £ g (n g ). 
Here the quadrupole quadratic forms A 2 ^ are given by 



A20 = A7 



C 22) 



.4 



22 



1 

Jf 



2x20 - v 2/7(j22 (2:20 + x 21 - 2x 
-1-^x22x42^ - (l/3\/77)<244 (I0x4 
2x 22 + \/277g22 (2x20x22 + VW2X21) 

+ (2/3Vl~4)<724 (x 2 2X40 + \/5X2lX41 + VT5X20X42 
+ v/2/33<744 ( 2x 42 X44 + 3X41X43 + 3^10/7x40X42 



(V3 X20X40 



10X21X41 



28x44 



35X21X43 



70X22X44) 



(5/V7> 



41 



(25) 



(26) 



5 



and the hexadecapole quadratic forms by 



A 40 = 



2x i0 + y / 2/35h 2 2 (3^20 — + x 2 22 ) — (2/\/77)/i24 [2\fhx2oXAn + 



• 1:5 ' 14x 44 ) 



1 



-SV2,x 2 2Xi2 J + ^2/1001/144 (9a;| + 9x 2 41 - llx 2 42 - 2lx 2 

2^42 + (l/V7)h 22 (V6x 2 ox 2 2 - 2x%^\ + y/3/385h 2i 
x ^6\/5a; 2 22:4o + 9x 2 ia;4i - (8/\/3)x 2 oX4 2 - 5V7x2ix i3 , + 2\f\ix 22 x^ 



Aa 4 = —r 



1 

St 



X42X44 + 2v / 7a;4i^43 — Hx^qx^ — Qx\ x 
2x A a + h 22 x\ 2 + (2/V55)/i 2 4 {VQx 22 Xi 2 + V21x 2 ix i3 + 2^X20^44 
+ (l/Vl43)/i44 ( 2V\AxioXa + 2v / 35a;4iX43 + 3^5^ 



The one-body expectation values in Eq. ( p5| ) are given by Eq. (reflbsd) and 

, N 



Finally, from Eq. (110) the cranking term is 



2^41 



2x^2 



2x 2 3 



2^44) 



2N r 
77 



£21 



(^\/6a; 2 o + 2x 2 2^j + £41 ^2\75x 4 o + 3V2x 42 ^ + x i3 , (^/Tlx^ + 2\/2x^ 



(27) 



(28) 



(29) 



In view of the numerous parameters in the sdg-IBM, we limit our discussion to a realistic range tailored to data 
p3| . Accordingly, the quadrupole parameters {^22, 924, 944} are scaled from their SU(3) values with a single factor 
q as suggested by microscopies [ JL4| . The hexadecapole parameters {/122, ^24, ^44} are determined from those of qji 
through the commutation condition 



[h, q] = 0, q jt = (j0l0\20)q jh h jt = (jOlO\40)h ju 



(30) 



which reproduce the available E4 data reasonably well |l3). We adapt the realistic sd-IBM parameters used in Fig. [I], 
k 2 = —20 keV, q = 0.5 and td = 1.5Nk,2, and use further, e g = 4Nk 2 and K4/K2 — 0-0.4. Energy levels are very 
well described in the sdg-IBM , therefore we do not dwell on a study of MOI apart from noting that all the SCC 
calculations of exhibit the characteristic rise with increasing frequency seen in experimental spectra. We focus 
instead on the shape question which could not be addressed in the 1 /N expansion approach. With the parameter set 
described above, the SCC calculations of 7 in the sdg-IBM are very similar to those of the sd-lBM with the realistic 
set of parameters (Fig. right). Because the hexadecapole interaction is a relatively novel feature not used in the 
sd-lBM, we briefly comment on its effect. Inclusion of the hexadecapole interaction does not cause any deviation 
from the axial shape but depending on its sign, it cither decreases the fluctuations in 7 (attractive) or increases them 
(repulsive) by a few degrees. This is very similar to the effect of changing the strength of the quadrupole interaction, 
in that, a larger k 2 leads to a deeper energy well in the 7 direction, therefore reduces the fluctuations, and vice versa. 
In certain parametrizations where the boson energies are neglected, it is possible to obtain non-axial shapes in the 
sdg-IBM |15| . However, once a g-boson energy of e g ~ 1 MeV is included, as demanded by the experimental spectra, 
such deviations from axial symmetry quickly disappear. In short, the congruence found between the sd- and sdg-IBM 
results with regard to the shape variable 7 is a direct consequence of the fairly high g-boson energy, which limits the 
effect of g bosons to a perturbative range at low to medium spins. 



V. CONCLUSIONS 



We have tested the self-consistent cranking method using the exact IBM results obtained from dynamical symmetries 
and numerical diagonalization. The SCC results for the dynamic MOI are in good agreement with the exact ones, 
especially in realistic cases where the deviation is at most a few percent. The study of shapes using SCC, on the other 
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hand, remains problematic due to the generation of spurious frequency-dependent terms in shape variables, which 
become dominant at high frequencies. Nevertheless, SCC can give reliable information on evolution of shapes in the 
low to mid frequency range. Our study of shapes using SCC indicates that both the sd- and sdg-lBM with one- and 
two-body interactions retain the axial shape (7 = 0) with fluctuations around 10°, as long as the model parameters 
are restricted to a realistic range that reproduce the data. The only frequency dependence in shape variables is 
observed in the fluctuations in 7 which is slightly reduced with increasing frequency. These results suggest that if the 
triaxial effects are important and need to be included in the IBM, then the remedy should be sought in three-body 
interactions rather than g bosons. 



APPENDIX A: THREE-BODY OPERATORS 

Here we derive the expectation value of the three-body quadrupole interaction in Eq. (|l4|). The scalar product of 
the three quadrupole operators is defined as 

Q ■ Q ■ Q = [QQ] (2) Q. (Al) 
The expectation value of this three-body term in the condensate state Eq. (||) is given by 

(Q ■ Q • Q) = Yl (-^n^M^)(0\b N Q^Q^Q-^) N \0). (A2) 

fJ.lfJ.2l-l 

Writing the quadrupole operators explicitly and commuting all the boson creation operators to the left, one obtains a 
three-body term, 3 two-body terms and a one-body term. Using boson calculus, the matrix elements of these normal 
ordered operators can be evaluated in a straightforward manner with the result 

(Q-Q-Q) = ]T (-l)"(2 Ml 2 M2 |2/x) Yl (-l) ni+n2+n3 Qhh<lhh<lhh 

/J1/J2A 1 31J2J3 ! 1 ! 2 ! 3 

mim 2 m 3 n 1 n 2 n 3 

x(jimilin 1 \2^ 1 )(j2m 2 l2n2\2[J.2)(j3rn 3 l 3 ri 3 \2 - /i) 

x { IV (IV 1)(A^ 2 l )xj irni Xj 2rn2 Xj 3m . 3 xi 1 — ni xi 2 — n2 xi 3 — n3 

-\-IVyIV \')\Xj imi Xj 2rrl2 Xi 2 — n2 Xi 3 — n3 8j 3 i l 5 ra3 — ni -\- Xj l7ni Xj 2m2 

XX[ 1 ~ ni Xi 3 — rL3 5j 3 i 2 5 rn3 — n2 -(- Xj 1 mi ^j^ms •El2—7i2^ls—n3 ^72^1 ^2 — ^1 ] 
-j-NXj l7ni Xl 3 — rl3 Sj 2 i 1 S m2 — ni Sj 3 l 2 S rrl3 — n 2 }. (A3) 

The three-body part (the first term) is seen to involve a triple product of the one-body expectation values A 2/J defined 
in Eq. (||). In the effective two-body terms, three C-G coefficients can be summed to yield a 6 — j symbol and a C-G 
coefficient. After some manipulations of the summation indices, all 3 terms can be shown to be equivalent. Finally, 



the sum of the four C-G coefficients in the effective one-body term gives a 6 — j symbol. Thus, Eq. (A3) can be 
written compactly as 

(Q-Q-Q) = N(N - l)(N - 2) £ (-l)^(2 Ml 2^ 2 |2^)A 2w A 2M2 ^ 2 _ M 

+3A^ - 1) ^(-l)^i 2M A 2 _ M + N^^L, (A4) 



where we have introduced 

i V. J 

25 ^ f / i 2 1 

£l = 2TTT ^ 1 2 2 j (A5) 

ij 

and A 2a , is obtained from A 2/J by replacing qji — > qji in Eq. (0). Note t hat symmetry of qji is retained, i.e., qji = qij 



Using the fact that A 2i = A 2 -i = and A 22 = A 2 _ 2 , the sums in Eq. ( |A4| ) involving A 2/i can be carried out to yield 
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(Q-Q-Q) = N(N- l)(N - 2)^/2/7(-A 2 20 + 6A 2 22 )A 2Q 



+3N(N - 1) (i 20 ^2o + 2A 22 A 22 + JVj^ ej: 



lm , 



(A6) 



lm 



In the case of the sd-IBM, A20 and A22 are given in Eq. (|T^). The contracted parameters in Eq. (Ac) become 

922 = 1 - 3 X 2 /14, 



902 = X, 
eo = 5 X , £2 = 2x - 3 X 3 /14- 



(A7) 



Using these values of qr^ , one obtains for A20 and A22 



A20 = A7 
i22 = A7 



2 X ^ - V / 2?7(l " 3x 2 /14)(x 2 + x\ - 2x\) 
2 X x 2 + ^277(1 - 3 X 2 /U)(2x x 2 + ^3/2x1) 



(A8) 



In the sdg-lBM, these parameters are given by 
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14 9 22 + y9 2 4 ; 924 
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£0 = 5922, £2 = 2(J 2 
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FIG. 1. Comparison of the cranking calculations of the dynamical moment of inertia in the sd-IBM (lines) with the exact 
results obtained from the dynamical symmetries SU(3) and 0(6) (filled circles), and numerical diagonalization (x = —\/7/A 
with e = and e = 1.5Nk) (open circles). In all cases, N — 10 and k = —20 keV are used. 



FIG. 2. Boson number dependence of the cranking dynamical moment of inertia in the SU(3) and Q(6) limits of the sd-IBM 
(circles). The lines that trace the circles are obtained from Eq. (^lj) in the SU(3) case and from Eq. ( [22] ) in the 0(6) case. 
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FIG. 3. Cranking calculations of the triaxiality angle 7 in the sd-IBM for the three cases shown in Fig. 1 for N = 10, 15, 20. 
The solid lines indicate average 7 obtained from Eq. (^) and the dashed lines correspond to its fluctuations given by Eq. (|l4|). 
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